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Abstract 

We prove the Central Limit Theorem for finite-dimensional vectors of linear eigenvalue statis¬ 
tics of submatrices of Wigner random matrices under the assumption that test functions are 
sufficiently smooth. We connect the asymptotic covariance to a family of correlated Gaussian 
Free Fields. 


1 Introduction 

The goal of this paper is to prove the central limit theorem for the joint distribution of linear 
eigenvalue statistics for submatrices of Wigner random matrices. 

Let and {Wjk}i<j<:k<n be two independent families of independent and identically 

distributed real-valued random variables satisfying: 

E[Wjk] = 0, E[lCfc] = 1 for j < k, and = T. (1.1) 

Set W = {Wjk)'j k=i with Wjk = W^j- The Wigner Ensemble of normalized real symmetric 
n X n matrices consists of matrices M of the form 

M = -^W. (1.2) 

yn 

For a real symmetric (Hermitian) matrix M of order n, its empirical distribution of the eigen¬ 
values is defined as jiM = ^ ELi where Ai < ... < An are the (ordered) eigenvalues of M. The 
Wigner semicircle law states that for any bounded continuous test function (p : M —>• R, the linear 
statistic 

n 

- X = -Tr(</^(AT)) =: trn(vj(M)) (1.3) 

n n 

2=1 

converges to f ip{x)dnsc{dx) in probability, where Hsc is determined by its density 

^(^) = ^T4-x21[_2,2](x), (1.4) 

see e.g. [I3], [5], [I]. 
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2 Statement of Main Results 


For a generic random variable in what follows denote by ^ For a finite set 

B C {1,2, ...,n} denote by M[B) the submatrix of M formed by the entries corresponding to 
intersections of rows and columns of M marked by the indices in B, which inherits the ordering. 
For example, 

=(«::«::) ^ 

Let Bi, • • • , be infinite subsets of N such that Bi, 1 <i < d, and their pairwise intersections 
have positive densities. Denote 

B'^ = Bin{l,2,... ,n}, l<i<d, 

Hi = \B^\, 1 < i < d, 
nim = \Bf ^ B^ 1, 1 < f < m < d. 

We assume that the following limits exist: 

7 ; := lim — > 0, ^im '■= lim 1 < I < m < d. (2.5) 

n^oo Ti n—)-cx> fi 

If it does not lead to ambiguity, we will omit the superindex n in the notation for 1 < i < d. 

I D| 

For an n X n matrix M and B C {1,2,... ,n}, consider a spectral linear statistic Xl=i 
where {Az}J^{ are the eigenvalues of the submatrix M{B). We are going to study the joint fluctua¬ 
tions of linear statistics of the eigenvalues. It will be beneficial later to view the submatrices from 
a different perspective. Consider the matrix = diag{P^), which projects onto the subspace 
corresponding to indices in B, i.e. 

= <n. (2.6) 

Define 

mB pB^pB^ ^2.7) 

n 

= ( 2 - 8 ) 

i=i 

where {AP})L]^ are the eigenvalues of . Note that the spectra of and M{B) differ only 
by a zero eigenvalue of multiplicity n — \B\. As a result, when we consider the linear statistics of 
their eigenvalues the extra terms {n — |i?|)<^(0) cancel when we center these random variables. In 
general, when considering multiple sequences Bi, to simplify the notation write 

Mil) M^\ PW := P^\ A/-«[(^] := . (2.9) 

Also, denote by pdP the matrix which projects onto the subspace corresponding to the indices in 
the intersection Bid] B^., i.e. 

pdP = pil)p{r) ^ ( 2 . 10 ) 

Recall that a test function ip belongs to the Sobolev space Ps if 

/ CXD 

ii + \t\r\m\^dt < oo. (2.11) 

-CX) 

First we consider Gaussian Wigner matrices. 


( 2 . 2 ) 

(2.3) 

(2.4) 
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Theorem 2.1. Let W = {Wjk '■ Wjk = Wkj}^).^^ be an nx n real symmetric random matrix with 
Gaussian entries satisfying (dip and M = n ^1‘^W. Let Bi,.. .B^ be infinite subsets o/N satisfying 
/i2.3^2.5\} . Let ifi, ■ ■ ■ , (/jrf : M —)• M 6e test functions that satisfy the regularity condition ||(/Jz||s < oo, 
for some s > |. Then the random vector 


( 2 . 12 ) 

converges in distribution to the zero mean Gaussian vector {Gi, • • • , Gd) G with covariance given 
by 


Co^r{Gl,Gp) 

2 ^ oo X 

4 2 ^ ^ WtHi 

2 
TT 


/ / *’!(" + 7 ) ("’+») S 


lip 


— zw 


+ 


kP = ll \w\'^ = Ip 
Jm z > 0 Jm w > 0 

'yip{a^-2) XtiW 


lip - zw 


L 


d\ I 


^T^'^lilp J-2^ a/ 4'yi - A2 J-2^ a/47p - ^ 



dzdw 


(2.13) 


In the expression for the covariance, {g^i)k denotes the coefficients in the expansion of <pi in the 
(rescaled) Chebyshev basis, i.e. 


and 


Ti{x) = ^(y5z)fcT^'(x), {ipi)k 
k=0 


2 

TT 





dt 

\/•ill - 


(2.14) 


T2 (x) = cos 


arccos 



(2.15) 


Note the form of the kernel in the above contour integral expression for the covariance. Since it 
is the Greens function for the Laplacian on El with Dirichlet boundary conditions (appropriately 
scaled), we note that the limiting distributions form a family of correlated Gaussian free fields. This 
is consistent with the previous work of A. Borodin in [S], m for the covariance of linear eigenvalue 
statistics corresponding to polynomial test functions. 

Now we formulate our result for the non-Gaussian Wigner matrices. 


Theorem 2.2. Let W = {Wjk)l^,^i be an n x n random matrix and M = . Let Bi, ■ ■ - Bd 

be infinite subsets o/N satisfying i2.^2.4\ ) and I12.5\) . Assume the following conditions: 

(1) All the entries ofW are independent random variables. 

(2) The fourth moment of the non-zero off-diagonal entries does not depend on n: 


= K{W%] 

(3) There exists a constant cjg such that for any j, k, < erg. 

Let (/?!, • • • , : M —)• M be test functions that satisfy the regularity condition < 00 , for some 
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s > 5.5. Then the random veetor 112.1S\) converges in distribution to the zero mean Gaussian vector 
(Gi, • • • , Gd) S with covariance given by 


Cov(G^Gp) = Cov(Gz,Gp) + 


27,-A 2 /-W 27P-//2 

2'7r^7, 7p 7-2 v7F 47, - J-2^ V 47p - T' 


-.dfi 


where Cov{Gi,Gp) is given by 112.1,^) . 


(2.16) 


In the course of the proof of theorem 12.11 it was necessary to understand the following bilinear 
form. 

Definition 2.3. Let M be a Wigner matrix satisfying 0, and let be the projection 

matrices defined in (j2.6p and (j2.10j) . For functions /, 5 G Ps-,s > define 

{f,g)lr := lim i V E [/(M«),fc • 

n—>-cxD n ^^ 


j,keBinBr 


= lim —E 

n^oo Ti 


Tr 




(2.17) 


Remark 2.4. The bilinear form (•, •),r is well defined on %s x TLs as a consequence of proposition 
The bilinear form is also well defined for polynomial / and g, see section 13.21 and also lemma 
12.51 below. 

The following diagonalization lemma is an important technical tool for the proof of theorem 

[2T1 


Lemma 2.5. The two families of resealed Ghebyshev polynomials of the 

second kind diagonalize the bilinear form HS.lSl]) . More precisely, 


VWh 




kq 


/ \ ^+1 

\V^ 


Let f,g^Ps, for some s > ^. A consequence of 112.18]) is that 


(2.18) 


1 GyGl /■ 2 yTV 

{f,a)ir = Y^ - / / f{x)g{y) 

J-2^i J-2^ 




.k=0 


k/2 fc/2 
7 ,' Ir 


\/ 47 , - x^\/ 47 ^ - y'^dydx. 

(2.19) 


In sectionit will also be proved that, with f,g given as above, almost surely 


lim ^Tl' • 5 (mW)pW} 

00 


1 




rW'yi f 

J-2./^ J- 


4vr27,7r 7-2y=7 J-2./^ 


f{x)g{y) 


,k+l 


,fc =0 7 , Ir 


\J 47 , - X^\/Mfr - y^dydx. 

( 2 . 20 ) 
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Remark 2.6. Recall that the rescaled Chebyshev polynomials of the second kind are orthonormal 
with respect to the Wigner semicircle law, i.e. 


1 - 

/ U//{x)U/l{x)\/4'y x‘^dx-6kq. 

271-7 J-2^ 

(2.21) 

C/72v/7cos(0)) = 

sm(tl) 

(2.22) 


The paper is structured as follows. The proof of theorem 12.II appears in section 3 and the proof 
of theorem 12.21 appears in section 4. 


3 Proof of Theorem 2.1 

3.1 Stein-Tikhomirov Method 

We follow the approach used by A. Lytova and L. Pastur in [TO] for the full Wigner matrix case. 
Essentially, it is a modification of the Stein-Tikhomirov method. This approach was also used to 
prove the CLT for linear eigenvalue statistics of band random matrices in |9] , which is connected 
to our work through the Chu-Vandermonde identity. See section 13.21 While several steps of our 
proof are similar to the ones in m, the fact that we are dealing with submatrices introduces new 
technical difficulties. 

We will prove Theorem [2T] in the present section and extend to non-Gaussian Wigner matrices 
later. The following inequalities will be used often. A consequence of the Poincare inequality is 
that, for differentiable test functions (p, 

VariTryjfM)} < LUe \Tr{ip'(M)(ip'(M))*}] (3.1) 

< 4((j 2-M) ('sup|v9'(x)A . (3.2) 

\31G]R / 

See [To] for a reference. The next inequality is due to M. Shcherbina, see [H]. Let s > 3/2 and 
ip G Rs- Then there is a constant Cg > 0, so that 

Var{Trc/,(A/)}<C',||(/.||^ (3.3) 

Let e > 0 and set s = | -|- e. Recall that the regularity assumption on the test functions is that 
||(/?;|| 5 / 2 +e < oo, for 1 < / < d. There exists a > 0 so that 

Var{AA«[(^z]} = Var{Tr(^KM(R/))} < CM\l/2+e- (3-4) 

The inequality holds because of (|3.3I) . noting that M{Bi) is an ordinary \Bi\ x \Bi\ Gaussian 
Wigner matrix. Note that this bound is n-independent. 

R is sufficient to prove the GLT for all linear combinations of the components of (|2.12p . Gonsider 
a linear combination ^ := aiAf^^'^°[pi], and denote the characteristic function by 
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Z„(x)=E[e*"«]. 


( 3 . 5 ) 


It is a basic fact that the characteristic function of the Gaussian distribution with variance V 
is given by 


Z(x} 


g-^V/ 2 ^ 


(3.6) 


As a consequence of the Levy Continuity theorem, to prove theorem 12.11 it will be sufficient to 
demonstrate that for each x G M, 


lim Zn{x) = Z{x), 

n^oo 

where Z{x) is given as above with 


V := lim 

n-^-cx) 


d 

j;a?Var(AA«°[^,])+ 2 

1=1 l<l<r<d 


(3.7) 


(3.8) 


So V is the limiting variance of It will be demonstrated that Zn{x) converges uniformly to the 
solution of the following equation 


Z{x) = I-V [ yZ{y)dy. (3.9) 

JD 

Note that (13.6h is the unique solution of (13.9h within the class of bounded and continuous 
functions. Therefore, to prove the theorem, it is sufficient to demonstrate that the pointwise limit 
of Zn{x) is a continuous and bounded function which satisfies equation (j3.9p . with V given by (j3.8p . 
Observe that 

d 

Z'^ix) = iE[Ce“«] = i (3.10) 

i=i 

Now it follows by the Cauchy-Schwartz inequality and (13.4p that 

d _ d 

\Zl,{x)\ < ^ Iaziy Var{A7(0[<^,]} < Const^ Iazi IIVJ;II5/2+e. (3.11) 

1=1 1=1 

Since Zn{0) = 1, we have by the fundamental theorem of calculus that 

Zn{x) = 1 + / Z'^{y)dy. (3.12) 

Jo 

Then to prove the CLT it is sufficient to show that any uniformly converging subsequences 
{Znm} and satisfy 

lim Zn^{x) = Z{x), (3.13) 

rim^oo 

and 

Z'n^x) = -xVZ{x). (3.14) 

A pre-compactness argument based on the Arzela-Ascoli theorem will be developed below, which 
ensures that the subsequences converge uniformly, implying that the limit is a continuous function. 
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The estimate \Zn{x)\ < 1, for all n, shows that the sequence is uniformly bounded. Generally we 
will abuse the subsequence notation by writing {n} for a uniformly converging subsequence. Since 
(13.lip combined with ||</ 3 z|| 5 / 2 +e < oo justify an application of the dominated convergence theorem 
in ()3.12p . it follows from (j3.13p and ()3.14p that the limit of Zn{x) satisfies equation ()3.9p . Therefore 
the pointwise limit (13.71) holds. We turn our attention to the pre-compactness argument, and will 
arerue later that (I3.13P and (I3.14p hold. Similar notation is used as in m- Denote by 

Djk := d/dMjk] (3.15) 

[/(0(t) ;= := (3.16) 

u«(t) := Tr{P«[/W(t)P«}, u«°(t) := u«(t) - E{u«(t)}. (3.17) 

For the benefit of the reader, what is needed is recorded below. Recall that is a unitary 

matrix, and writing Pjk ■= (1 + Sjk)~^, we have 

n 

\uji\<l, J]|17«|2 = 1 , ||17«|| = 1. (3.18) 

k=l 

Moreover, 

Ojl.C'iht) = ViilUMB,) (<’ • <(«) + f/S • <(()) . 

where ^ 

f*g{t) := f f{y)g{t -y)dy. 

Jo 

Applying the Fourier inversion formula 



poo 

V’lW = / e"^^(pi{t)dt, 

J — OO 

(3.21) 

it follows that 



Now dehne 

J —OO 

(3.22) 


en{x) := 

(3.23) 

Usine the Fourier reoresentation of the linear eigenvalue statistics in (13.101). 

it follows that 


d pQQ 

Znix) = i'^oii (pi{t)Y^^\x,t)dt, 

1=1 

(3.24) 

where 




Yji^\x,t) :=E {t)en{x) . 

(3.25) 

The limit of Yn^ {x 

t) is determined later in the proof. Since 



Yj}\x,t) = Y^\-x,-t), 

(3.26) 


(3.19) 

(3.20) 
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we need only consider t > 0. It will now be demonstrated that each sequence is bounded 

and equicontinuous on compact subsets of {x G M, t > 0}, and that every uniformly converging 
subsequence has the same limit implying (I3.13P and (13.141) . See proposition 13.11 

Let ^{x) = and note that sup^gjj \ip'{x)\ = \t\. Applying the inequality (13.2p to the linear 
eigenvalue statistic [tp ], we obtain 


Var{uW(t)} = Var{AA(')[(^]} < 4(a^ + 


(3.27) 


Now set ip(x) = and notice that 


d 




Using the inequality p3.ip and the fact that n ^ETr(M(*))^ < cr^ + 1, it follows that 

Tr{(^'(MW)((^'(MW))*} 




< ihhtilE 


n 


Tr{l + t2(MW)2} 


(3.28) 


< 4{a^ + l)[l + {a^ + l)t% 

Using the Cauchy-Schwartz inequality, the bound |e„(x)| < 1, (|3.27l) and (|3.28p . we obtain 

y«(x,t) < VarV2{nW(i)} < 2(^2 + l)y^\t\, (3.29) 


and also 






(3.30) 


Observe that 


A 

dx 


e-nix) = ienix)'^ar 


r=l 


Using the above derivative with the Cauchy-Schwartz inequality, (I3.4I) and (I3.27p . we have that 

d 


A 

dx 




= ij;a,E[u«°(t)A/-W°[^,]e„(x)] 

r=l 

d 

< Var‘/2{uS)(t)}^|a,| Var^/2|_;y^(r)[^^j| 

r=l 

d 

< Const ■ i*iE Wr\ ||<7’r|l5/2-|-e- 


(3.31) 


r=l 


It follows from p3.29p . the mean value theorem combined with p3.30p and (I3.31jl . and | |</?r-|l 5 / 2 -i-e < 

oo, that each sequence Yn\x,t) is bounded and equicontinuous on compact subsets of The 
following proposition justifies this restriction. 
































Proposition 3.1. In order to prove the functions Yn\x, t) converge uniformly to appropriate limits 
so that id. 24 \ ) implies it is sufficient to prove the convergence of Yn\x,t) on arbitrary 

compact subsets o/ {x S M, t > 0}. 

Proof. Let 5 > 0. Recall that the regularity assumption on the test functions (pi are 

f {l + \h\f~^^ipi{h)\‘^dh<oo, 

Jr 

i.e. that ipi € Us, with s = 5/2 + e. Using the Cauchy-Schwartz inequality, it follows that 

/(I + \h\)\ipi{h)\dh 
Jr 


< 


dh 


(l + |h|)3+^ 


/ / ii + \h\)^+^\m)\^dh, 

Jr 


which implies that 


\h\ ■ \(pi{h)\dh < oo. 


(3.32) 


(3.33) 


A consequence of the finiteness of the integral in (I3.33D . for each 1 < / < d, is that there exists 
a T > 0 so that 


2 (cr^ + 1)^/2 ^ loil / \t\ ■ \ipi{t)\dt < 6. 

J\t\>T 


(3.34) 


1=1 


Using ()3.24l) . we can write 

d 


Zn{x) = i [ lpi(t)Yd\x,t)dt + i V'oz [ ipi{t)Yd\x,t)dt. (3.35) 

z=i z=i d\t\>T 


Then (|3.35p . (I3.29p . (I3.34p imply that 



d 

i ^ai 
1=1 



(pi{t)Yd\x,t)dt 


< 


[ \m\ 


1=1 

J\t\>T 



d 

< 

2{a^ + 




1=1 

< 

6. 



• \yn\x,t)\dt 



|t| • \ipi{t)\dt 


(3.36) 


Notice that the estimate (|3.36l) is n-independent, so that in particular the estimate holds in the 
limit n —?> oo. Since 5 was arbitrary, this completes the proof of the proposition. □ 

This completes the pre-compactness argument, which allows us to pass to the limit in (j3.24p 
and in (|3.12l) . and conclude that Zn{x) converges pointwise to the unique solution of equation (|3.9p 
belonging to (^^(M), implying (13.7h . and hence the conclusion of the theorem. Now we show the 
limiting behavior of the sequences Yn\x,t) imply (I3.13p and (13.141) . Consider the identity 

+ i f Md)^rhMi^^dh. 

Jo 
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Apply this identity, noting that = 0, if j, k ^ Bi, to obtain that 




rt n 


= I 




° i,fe=l 


(3.37) 


Recalling that Yn\x,t) = IE Un^°{t)en{x) 


, and applying the decoupling formula for Gaussian 


random variables, it follows from (13.371) that 

Yjl\x,t) = i J] E[M«t/«(ti)e°(a;)]dti 
2i /■* 

n Jo 


+- 


icr 


Y1 ^{xkeBi} IE 

l<_7</c<n. 

2 rt ^ 


D,kUfi{h)el{x) 


dti. 


n 


El 

i=i 


ii&Bi} 


E 


D,,uf]{B)el{x) 


dU. 


(3.38) 


It will be useful to rewrite (13.381) as 


rt n- 


n In ' 

•^0 j,k=l 


y«(x,t) = Y.hjMBt}a + Sjk)^ DjkU^iiti)e°n{x) 


dti 


+ 


i{a‘^ — 2 ) 


n 


rt ri 

Is 


=:Ti 


l{i 6 i?i}lE 


DnUfJ{t^)el{x) 


dU 


= .T2 


(3.39) 


The reason for the rewrite is that it splits the functions Yn\x,t) into a part that depends on 
the distribution of the diagonal entries and a part that corresponds to the same term as for the 
Gaussian Orthogonal Ensemble, for which = 2. Recalling that en{x) is given by (|3.23p . again 
writing /3jk = (1 + 5jfc)~^ and using the identity 

D,kTi-f{M) = 2f5jkf{M)jk, 


it follows by a direct calculation that 

d 

Djkenix) = 2if3jkxen{x) ^ ■ (3-40) 

r=l 
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Then for 1 < I < d, using (|3.40l) and (|3.19l) . it follows that 

rt rt\ 


Ti = 


-If 

n Jo Jo 

1 /■* /■*! 
nJo Jo 

r 

n Jo 


E 


E 


-t2)e°(x) 

_j,k=l 
n 

Td)u \TTd)l 


dt2dti 


dt2dti 


Y Mj,keBt}U}J{t2)U^J{ti - t2)e°^{x) 

j,k=l 

d 


j,k=l 


r=l 


jk 


dti, 


(3.41) 


and also that 


To = 


-( ct 2-2) /■' /■'! 


B Jo Jo 


If 


E 


Y ^{j&Bi}Uf {t 2 )uf {h - t 2 )el{x) 
i=i 


dt2dti 


= :T2i 


(^2 - 2)x 


n 


f 


E 


n a 

j=l r=l 


dti . 


= '-T22 


Using the semigroup property 

it follows form (j3.4ip that Ti can be written 

1 /“* 


(3.42) 


Ti = — 


n 


E 


0 JO 


ullfh - t2)u!I\t2)e°^{x) 


dt2dti 


= :Tii 


n 


tiE 


'0 


ulfiti)e°^{x) 


dti 


2x 


= .Ti 2 
d. rt 


/ E Tr{P^^^U^^\ti)P^^’^^(pl{M^^fP^^'^}en{x) 


B -! JO 

r=l 


dti ■ 


(3.43) 


= :Ti3 


Define 


S<'>(i) := i|E|««(*)|. 


(3.44) 


The following proposition presents the functions Yffx, t) in a form that is amenable to asymptotic 
analysis. 
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Proposition 3.2. The equation Yn\x,t) = Ti + T 2 , can be written as 

+ 2 /" [ vHHh - t 2 )Y^^\x,t 2 )dt 2 dti = xZnix) A^^\t) + Q®(t) + r^\x,t), 
Jo Jo L J 


where 


A^^\t) :=-2'Var / -E 
'0 L 


r=l 








n 


'^Otr j 

r=l ,=1 


E 


(3.45) 

dti, (3.46) 

U^’^{tMiM^%,}\dh, (3.47) 


and 


r^n\x,t) = 


— [ tiYj;’'\x,ti)dti 
n Jo 


1 /■*! 

Jo Jo 

^ r 

r=l Jo 


4o 40 

n 


u'J>^h-t2)u^>{t2)eUx) 


dt2dti 


rt rti 


E 


^ ./o JO 

x{a^ — 2 ) 


-2) r r 

Jo Jo 

d pt 


Elb-6sat^i?(^2)C/j?(ti-i2)e°(x) 

i=i 


(3.48) 

(3.49) 

dh (3.50) 

dt2dti (3.51) 


^j&BinBr}^ 


r=l j=l 

Proof. Begin with the term Tn, defined in ()3.43l) . Write 
1 7* 7*1 


U^'^{t,U{M^%jeUx)]dh. (3.52) 

(3.53) 


Tn = — / E 

nJo Jo 


u^J^°{ti - t 2 ) + nvn{ti -t 2 )) ■ ( uJ^°{t 2 ) + nVn{t 2 ) ) 6 ° (x) dt 2 dti, (3.54) 


so that 


Til = 


t ptl 



E 


'n Jo Jo 

rt rti 


{ti - t2)u^n^° {t2)el{x) 


dt2dti 


[ [ Xn(tl -i2)E L®°(t2)e°(x) 

Jo Jo 1 

[ [ hn(t2)E u^^)°(ti - t2)e°(x) 

Jo Jo L 

I ~ E [en{x)] dt2dti. 


dt2dti 

dt2dti 


=0 


(3.55) 
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Noting that 


E 


«n^°(^ 2 )e°(a:)j = Y^^\x,t2), E \u!p°{ti - t 2 )e°(x)j = - t2), 


and also that 


t pt\ 


[ [ Vn{t2)Yjl\x,ti-t2)dt2dti= j I Vn{tl -t2)Y^''\x,t2)dt2dti, 

Jo Jo Jo Jo 


it follows that 


Til = 


1 /■* /■*! 



E 


./o ^0 

rt pti 

-2 

lo Jo 


ul^>°{h - t2)ul>°it2)e°^ix) 

[ [ Vnih-t2)Y^''Hx^'>^2)dt2dti. 

Jo Jo 


dt2dti 


(3.56) 

(3.57) 

(3.58) 


The term ()3.56j) goes into the remainder, which becomes (j3.49p . Also, (I3.57jl is added to the 
left-hand side of (j3.45|) . Now consider the term T 12 , defined in (I3.43|) . We have that 


Ti 2 = - 


1 

-- / tiYji^\x,ti)dti, 

n Jo 


(3.59) 


which becomes p.48l) in the remainder. Consider the term T 13 , also defined in ()3.43p . Writing 
J rt 


Ti3 = - 


2t 7* r 1 

— Yar E (enW + ^n(x)) dh, (3.60) 

^ r=l do L J 


it follows, with An\t) given by (13.461) . that 
Ti3 = 

^ /*t r 

E YT{pd'^Ud\ti)pd^^'^ip'^{MYPYe°n{x) 

“1 JO L 


2x 

n 


dti 


r=l 


-h xZn{x)Ad){t). 


(3.61) 

(3.62) 

(3.63) 


Then (I3.6ip becomes (j3.50l) in the remainder, while (13.621) remains on the right-hand side of 
(1051) . 

Now consider the term T 21 , defined in (j3.42p . This term becomes (j3.5ip in the remainder. 
Finally, consider the term T 22 , also defined in (j3.42p . Write 


T 22 = — 


n Jo 


E 


• ielix) + Znix))Yotr V;(M^^^)P« 

j=l r=l 


n 


dti, 

(3.64) 
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so that, with Qn\t) given by (|3.47p 


T 22 = 


(fj^ — 2)x 


n 


/ f n d 

j=l r=l 


33 


+ xZn{x) ■ Q^n\t)- 


dti (3.65) 

(3.66) 

(3.67) 

The term (I3.65P becomes p3.52p in the remainder. Also, the term (I3.66P remains on the right-hand 
side of P3.45I) . This completes the argument for proposition 13.21 □ 

We now turn our attention to the remainder term, rn\x,t), of proposition 13.21 The content of 
the following proposition is that the remainder is negligible in the limit. 

Proposition 3.3. Each term of rn\x,t) converges to 0 uniformly on compact subsets of {x G 
M, f > 0}, for 1 < I < d. In other words, we have the uniform limit 

lim r^\x, t) = 0. 


(3.68) 


Proof. Begin with the term (j3.48p . Applying the estimate (I3.29|) . we obtain 


- f tiYj^’-\x,ti)dti 

n Jo 


< -r 
n 




^ 2 (u^ + l)V^ 3 


n 


= 01 - 

n 


(3.69) 

Now consider the term (13.491) . Using the bound |e°(a;)| < 2, the Cauchy-Schwartz inequality, and 
(I3.27P twice, it follows that 

X P 


n 



E 


0 JO 


- t2)u^n'’° {^2)el{x) 


dt2dti 


< -t 2 VarV 2 {^«(t)}Vari/ 2 { 40 (i)} 


n 


^ 8 (a^ + l)V ^ ^4 


n 


= o{- 

n 


(3.70) 

Consider the term ()3.5np next. Applying (j2.19p of lemma [ 2 A] to the exponential function and 
(/?(,, and noting that (/j), G it follows that 


lim —E 
n—)-oo 72 


Tr 


I p(0 [/(O (t 1 )pO:U (M W )p('')} 

1 r'^VTi 

A 2 — / / £’“>;(!/) 

in2^nr 3-2^ 


EU‘(i)crfe) 


Lfc=0 


kl2 k/2 
h Ir 


1 / 47 ; - x 2 - y'^dydx. 

(3.71) 


14 













































While the exponential function does not belong to "Hs , we can truncate the exponential 

2~'~ 

function in a smooth fashion outside the support of the semicircle law, so that the truncated 
exponential function belongs to "Ha.. We may replace the exponential function by its truncated 
version because the eigenvalues of the submatrices concentrate in the support of the semicircle law 
with overwhelming probability. Then 


lim -Tr IP(')(t 1 ) ip' )PW \ 

n—>oo n I j 

4 vr 27 ; 7 ,. J_2^ J-2^ 


itix J 


^'Ay) 


ET‘(i)Tfe) 


, fc =0 


k/2 k/2 


■sjA'ji - - y‘^dydx. 

(3.72) 


Here it is not so important to know the exact value of the limit, but we will use the fact that 
we have convergence in the mean and almost surely to the same limit. Note the convergence in 
(j3.7ip implies that the sequence of numbers 


1 . 


-E 


n 


Tr 


I pd) [/(O (i ^ W )pM I 


is bounded. Also the convergence in (j3.72p implies that the random variables 

1 , 


-Tr 


n 


I p(') (t 1 )p w I, 


are bounded with probability 1. Using (I3.71|l and (j3.72p with the dominated convergence theorem, 
it now follows that 


lim E 

n—^oo 




n n 

Combining the bound |en(x)| < 1 with ()3.73l) . it follows that 
1 


= 0. (3.73) 


-E 


< E 


Tr{P«[/(0(ti)p(7O(^(.(MM)pM}e°(x; 

-TrP(')pW(ti)p('’’'V(-(MW)pW - -E|TrP(')pW(ti)p('’''V;(MW)pW|^ e„(a:) 


(3.74) 


Then, using (j3.74p in the remainder term (j3.50p . it follows that 

d ft 


‘lx V—^ 


E 


r=l 


Tr{PW p(') (t 1 )p(^Py,; )PW }e° (x) 


dU 


0 as n —>■ oo. 


(3.75) 
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Consider (|3.51l) . which is the next term in the remainder. Observe that, again using the Cauchy- 
Schwartz inequality and the fact that |e„(x)| < 1, 


E 


n 


i=i 


= E 


< E 




j&Bi 


j&Bi 


f E - *2) - E - fe) 


S Var‘/^ .j i E - *2) 


J&Bi 


. 


j&Bi 


(3.76) 


For fixed j, p,q G Bi, using ()3.19p , 

= 2il3„ fu'»(t-h)uf;;(h)dh. 


'0 


(3.77) 


Usine (13.771). recalline that B^r, = (1 + ^pq) < 1, and the Cauchy-Schwartz inequality, it follows 
that 

Dp,U^’^{t)\" < m f \uf^{t - h)uf;i{h)\^dh. (3.78) 

'J 0 

Using (13.781) . the fact that |C/j^^(t)| < 1, and the inequality 2ab < + 6^, it follows that 

Dp,{uf]{t2)uf]{B-t2)}^ 


< 2\DpqU'fht2)\^ + 2\DpqU))>{tl - t2)\ 


(i). 


33 

8|t| {^f^^\uf^{t2-h)uf^{h)\‘^dh + '"\uf)(t,-t2-h)uf^^{htdh 


33 


< 


(3.79) 


Using the Poincare inequality, (I3.79p . adding more nonnegative terms, and using the property 
of the unitary matrices that 

n 

k=l 


(3.80) 
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it follows that 


Var 




< Y. ^ [(Afg) 
p<q 

8 (a2 + l)|t| " 


E 




< 


< 


n i 

p=lq=l L-'U 

8 ((t^ + l)N f 




n 

^2 


jp \ ^ / jq 

2 rti-t2 


dh -j- 


uSiti-t2-h)ul^ihf 


jq 


dh 


Ufi(t2 - h) 


dh -|- 


'0 


Uf^{h-t2-h) 


dh 


< 16K + i)I^L ^ 


n 


= O 


(3.81) 


Now, combining (j3.76jl with (j3.8ip . we have that 


E 




n 


i=i 


= O 


(3.82) 


and it follows that 


(CT^ - 2 ) 


n 


t rti 


E 


0 JO 


^l{,6B,}^j?(t2)C/j?(ti-t2)e°(x) 

i=i 


dt2dti 


= O 


(3.83) 


Now consider the final term of the remainder, given by ()3.52p . We apply the identity below 

/ CX) 

hprih)U^^\h)dh, (3.84) 

-OD 

which is a consequence of the matrix version of the Fourier inversion formula (j3.21jl . Using (j3.84p . 
the finiteness of the integral ()3.33p . the above estimate p3.82p . and the dominated convergence 
theorem, we have that 


x((T^ — 2 ) 


rt ri 


n 




r=l -^0 


{jeBinBr} 


E 


U«(ti)(y.;(MM)),-,e°(x) 


dti 


— “ 2 )a;r 


r=l 


rt roo 1 p 

/ / hMh)- E E ^j?(U)C/j;^(/i),ye°(x) 

^0 J—oo ^ *- 


dhdU 


j&Bl^Br 


0 . 


(3.85) 


Combining (I3.69p . ^3.701) . p3.75p . ^3.831) . p3.85p . and comparing to the remainder term p3.48l) . 
the proposition is proved. □ 
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The goal now is to pass to the limit in (13.451) . In what follows let {[/^(x)} denote the (rescaled) 
Chebyshev polynomials of the second kind on [—2y^, 2^], 


i=o \ J / \ 


2^7 


k-2j 


(3.86) 


Proposition 3.4. Let An\t) be given by Qn\t) given by {3.4'^ , and Vn{t) given by Ii3.44\ ) 

. Then the limits of An\t),Qn\t) and Vn{t) as n —>■ oo exist and 


m. 


^«(t) := lim ^W(t) 


1 “ Q, /■* r'^VT /■2v/tv _ _ 

^ 47 / - - y^Firix, y)dydxdti, 

11 7r Jo J-2./^ 7-2,/^ 


^^^7« ^ 7r Jo J- 2 ^J- 2 ^ 


where 


Fir{x,y) = ^Ul\x)U2f{y) 


h 


k+l 


k=0 


k/2 k/2 ■ 

h Ir 


(3.87) 

(3.88) 


the limit of Qn^ (t) is given by 
Q^^\t) := lim Q^J^it) 


(cj2 -2) 'jirar /■* 


^^^11 ^ 7 r Jo J-2^ 


(^(.(/r)-\/47r - ijfdyL, 

- 2 y> 

(3.89) 


and the limit of Vn{t), after rescaling by 'yi, is given by 

1 1 /* 2 y^ 

vJ\t) := — lim Vn{t) = - / e**^\/47; — x‘^dx. (3.90) 

7, n^oo 2'K-fi J-2^i 

Proof Recall that aJJ {t) = Jq [T\'{pW[/«(ti)P(^’^)v?;(MW)pW}] dh. In the full 

Wigner matrix case one has An{t) = —2 ^ETr{e^*^(p'(M)}dti, and the limiting behavior follows 

immediately from the Wigner semicircle law. In the case of submatrices with asymptotically regular 
intersections there are additional technical difficulties due to the fact that for the nxn submatrices 
Mil) ^ pii)MPd)^ we have 

Tr{pT'-)p(0(f)(^;(MM)p('’^)} = J] U^2it)Tr{M^%k, (3.91) 

j,keBinBr 

so that the summation is restricted to entries common to both submatrices, i.e. to j, k £ Bi D Br¬ 
it follows from lemma [23] that the limit of An\t) exists and equals 

7l«(t) = -2 Vo, )ir dti, (3.92) 

r=l -^0 
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where 


1 r‘^y/Ti /■2v/> ,_ ,_ 

= —^- / / e^^^^ip'^{y)Fir{x, y)\/^7i - x‘^V‘^7r - y'^dydx. (3.93) 

4vr^7z7,. J_2^J-2^ 


This establishes (I3.87j) . The proof of lemma 12.51 will be given in section 13.21 

We turn our attention to Qn\t). First it will be argued that the variance of the matrix entries 
converge to zero. Using the Poincare inequality, (j3.78p . (j3.80p . and proposition 13.11 it follows that 


Var ( 


< ^ E [(M«) 

4(a2 + l)|t 


E 


Dp,U]’it) 


< 


< 


n 


/*Cl 

P=1 0=1 -'U 


dt2 


4((T2 + l)|t 


p=i q= 

n 


n 


I r^i 

p=l “'O 


dt2 


< 




,2 11.^2 
n 

Note that in the course of the calculation (13.941) . we showed that 

2 " 


p<<i 


DpqUf]{tl] 


< 4t?. 


The Cauchy-Schwartz inequality implies 

I {I + t\)\(pr{ti)\dti < . 


dt^ 


R(l + tf)V2+e \ 

Since ||¥?r|| 5 / 2 -i-e < oo, we have the estimate 

/ CX) 

t\\(pr{tl)\dti < OO. 

-OO 

Using the Cauchy-Schwartz inequality and (I3.84p . it follows that 




33 


L 


— oo 
oo 


ti(pr{ti)DpqUf^ {ti)dti 


< 


ti\ipr{tl)\dti- j \(pr{ti)\- DpqUf^{ti) dti. (3.98) 


(3.94) 


(3.95) 


(1-hti)^/2+^|^r(ti)pdti. (3.96) 


(3.97) 
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Using the Poincare inequality, (I3.95p . (I3.98p . we obtain 


Var 


< 


E® 


E 


p<<i 
^2 




33 


< 


< 


+ 1 ) 


4(a^ + 1) 


n 


/ cx> _ poo 

tl\!pr{ti)\dti I |^,.(ti)|E 

/ oo 

-oo 


DpgUgih) 


dti 


Using (13.971) . (I3.99p . p3.94p . and the Cauchy-Schwartz inequality, we obtain 

Cov{[/j')(U),(^;(MW),,} < ^Jvar{U^^{t^)} • = O (n'l) . 


(3.99) 


(3.100) 


Using p3.100p it is justified to replace the expectation E[ujj\t)(/ 9 (.(M(^))jj] by the product E[Uj’j{t)]- 
e[(^;(mW),-,], when passing to the limit. We use proposition 2.1 of m, which guarantees that for 

lim E[f{M)jj] = [ f{x)df3scix). (3.101) 

n^oo 

In order to apply this asymptotic to the exponential function, which is smooth enough, we truncate 
the function in a smooth fashion outside the support of fisc- We are justified in replacing the 
exponential function by its truncated version because the eigenvalues of the submatrices concentrate 
in the support of the semicircle law, with overwhelming probability. It is for this same reason that 
we may assume is compactly supported. This function is not sufficiently smooth, but we can 
avoid this problem by a density argument using standard convolution, and then apply the bound 
(|3.3I) on the variance of linear eigenvalue statistics. 

Let r] E C')?°(M) satisfy Jj^r]{x)dx = 1, and consider the mollifiers r]y{x) := y~^r]{xy~^). Then 
p>'j. * tfy E (M), and using standard Fourier theory it can be shown that 


lim|l¥^r-¥^r*^ylli/2+e = 0- 


(3.102) 


It follows from (I3.100p and (|3.10ip that 


lim — 1 

r—^r<r, n / ^ 


n^OO Tl 


{jeB,nBr} 


E 


c/«(t)c/,;(MW), 


33 


1=1 


r- 2 ^ 


'Jlr 


- \^d\ 




^27r-/iJ_2^ J y27r-frJ_2^ 

Using (I3.103p . we pass to the limit in (13.471) . and obtain (13.891) . The limit of 

= -E[u«(t)] ^E[Tr{P«C/«(t)P«}], 

n \ni\ 


(f'^{fi)Y^i'yr - fJ-'^dfi . (3.103) 
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is given by (rescaled) Wigner semicircle law, as a consequence of the zero eigenvalues. Alternatively, 
it can be computed using the bilinear form in lemma [2Al with f{x) = e*** and g{x) = 1. To facilitate 
solving the integral equation (I3.105p . below, it will be useful to rescale by ji. We obtain 




'll 


1 


r2^ 


2TT'yi J-2^ 

which establishes ()3.90p . The proposition is proved. 


1 / 47 ^ — x'^dx , 


(3.104) 


□ 


Now using propositions 13.2113.3113.41 we pass to the limit rim —)• 00 in (I3.45p . and determine 
that the limit of every uniformly converging subsequence satishes the equation 


Y\x,t) + 2'yi f f - t 2 )Y^’‘\x,t 2 )dt 2 dti = xZ{x) (t) + , (3.105) 

Jo Jo L J 

where AJ\t) is given by (I3.87jl . is given by (I3.89p . and vJ\t) is given by (j3.90p . 

Now the argument will proceed by solving the integral equation p3.105p . We use a version of 
the technique used by L. Pastur and A. Lytova in m, to solve this equation. Define 


f{z) := - 4ji - z)/2ji, 


(3.106) 


which is the Stieltjes transform of the rescaled semicircle law, where z^ — Yg = z -h 0(l/z) 
as z —>■ 00 . A direct calculation shows that = /, where denotes the generalized Fourier 
transform of We obtain 

r-oo 


1 1*00 

vJ\z) := V4t - 

2v^^7^ Jo J-2^/^ 

- 1 / 47 ; — x'^dx 


'-2YTi 
r'^VY 1 


2 vr 7 z 

= f{z). 


2ATi 


x — z 


We check that 

Set 


2 ; + 2-fif{z) = v^ 2;2 _ 4 ^^ ^ 0 , 3mz^^. 




e'-^^dz 


1 r2YYi 


L z + 2jif{z) ^ J-2^ a/ A-fi - A 2 


/ 


(3.107) 

(3.108) 

(3.109) 


after replacing the integral over L by the integral over [— 27 /, 27 ;], and taking into account that 
y/z"^ - 47 / is on the upper and lower edges of the cut. Then the solution of (I3.105P 

is 


YJHx,t) = -xZ{x) + 

Jo “^1 


dti. 


(3.110) 
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Then, with Fir given by (I3.88p . 


[ T{t 

Jo dti 


1 


2TT^'yi 

xFir{x, y)dydxdXdti 


t 1-2^^ 1-2^ f2^r 


y^ar r rv^' rv^' r 

^ 7r Jo J- 2 ^iJ- 2 ^iJ- 


-2V> 




1 


E 


CXr' 


r2FF r2y/Yi r2FF - e*‘^] (Pr{y) 


2i7r37; ^ 7 ^ V-2V7i -'-2V7i J-2FF 


a/ 47 ; - 


Fir{x, y)dydxdX, 
(3.111) 


and 
(■t 


f T{t-ti)-^QJ\ti)dti 
Jo ®i 




X 


^ 7r 7o .1-2^ -'-2^71 a/ 47 ; - A2 
_ 

/ (/?(, (//) 47r - yfldy,dti 

J-2^ 


'yiria'^ - 2 ) 

d 

Oir 

pFFi 

pFFi 

' pitr} _ pitX' 

47r37;i 


'-2FYiJ 

-2FTi 

rj — X 


f'^ 

V 47z - A2 J-2^ 


tp'r{yL)\/A'yr - fjJdfi. 

(3.112) 


Using the regularity condition ||(^z|| 5 / 2 +e < 00 for 1 < / < d, (I3.11ip . (I3.112p . and the dominated 
convergence theorem to pass to limit in p3.24p yields 


Z'{x) 


= i 


/ oo 

(pi{t)Y^'‘\x,t)dt 

t=i 

xZ{x) ^^aiar r 

r=i T7r 7-00 J-2^i J-2^ J- 


27r^ 


2y^ f‘2yj^ Tgitai — 

I I ipi{t)^ 




(x- A) 


X a;V 47 ^ y^ y)dydxdXdt 


x/Aji - A^ 


7ir(<7^ — 2)xZ(x) ^ ^ f2VY f2^ 


47r3 


EE 


X 


r=l -2x/Tl d-2y/Yi 

f ifr{y)\/4:'yr - yjdy.dt. 

J-2,/^ 


47 ; - 7/2 


dr]dX 


r}-X \ - A2 

(3.113) 
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Applying the Fourier inversion formula (13.211) . it follows that 


Z'(x) = 


d d 


aittr ^i{x) - ipi{X)]ip'j.{y) 


xZ(x) s^s^aiar f 

2^^ ^ ^ 7ar J-2^ J-2^ J- 


(x - A) 


a/ 47 / - 


xFir{x, y)dydxdX 


—xZ{x) 


'yirjcr'^ - 2 ) 

47r3 




V’liv) - fiW 


- 2^1 L ^ J a/47/ - A2 


^ 47 , - r /2 


drjdX 


X 


/ y::'r{y)\J- y^dy. 

J- 2 ,/Vr 


'- 2 v^ 

We will use the fact that 

K(i-) - r’(A)] dA 


/: 


TT 


'- 2 ^ {x - y^47 - A 2 2^/7 

Expand the test function (/j; in the Chebyshev basis to obtain 


k>l. 


J!2. 9 rh 

^i{x) = Y^{^i)kT^\x), {ipi)k = - ^i{t)T2\t) 

^ J- 2 .FT 1 V 47z - i 


(3.114) 


(3.115) 


(3.116) 


/c=0 


Returning to the computation of Z'{x), using (13.114^ . (I3.115p . and (I3.116p . it follows that 
4^ ,=1 7z ' It J-^V^ J-^FTr 


xFir{x,y)dydx 

, o d d 00 


SVT^ 


_ _ 3/2 

1=1 r=l h Jr k=l 


/ ^7^ii(r?)\/47z - 72^7 

, J-2^i 


r'^VF - 

/ ip'r{y)V^Jr - y‘^dy. 

J—2^/^ 


(3.117) 


F2. r'^VF _ 

/ _ C^fcii(^)\/47z - y‘^dy = 2 ^ 

fc=i 


-2V7i V47 z - A 2 


1-2^/^ 

Using the orthogonality of the Chebyshev polynomials (I2.2ip . 

^k-ivi)\^n->r^'i-^yni 

Integrating by parts yields 

r'^7(M)vswu</d=r'^ 

J-27> 

so that 

7Zr(o'^-2) A(/?z(A) 


dA. 


(3.118) 


'- 2 U> a/ 47 ^ - y? 


dy, 


J- 2^1 a/47z - A2 7-27> A/47r - 






(3.119) 


(3.120) 
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Since 


(3.121) 




we expand (pr{y) in the Chebyshev basis to obtain 


(Priv) = X] 

V m=l 


(3.122) 


Recalling that Fi^ is given by (I3.88p . it follows that 


EE m{(pi)k{y^r)r 


k=l m=l 
oo oo oo 


fwyr ^ ^ ! - . - 

/ / ^k-\^^')^Z-\^y)'J- a^V47r - y^Fir{x,y)dydx 

J-2,/y7 J-2,/y7 


/‘ 2 V> 


fc=lm=l7=0 J-2^J-2^ 


^i +1 

X V 47z - a;2^47^ - y^dydx . 

7/ Tr 


4vr^7z7r ^ k{ipi)k{(pr)k 


hr 


k=l 


(fe-l)/2 (fc-l)/2) 

.7/ 7r 7 


(3.123) 


Using P3.123D . (13.1181) . ^3.1191) and (I3.120p . in (13.1171) . it follows that 


'rViF . . 

rylf \ XZj[X) 

Z [X) = - 


1=1 r=l 


(cr^ - 2 ) 'Jir 
2 


CX) / 

(7’z)l(7^r)l +'^k{lfl)k{‘Pr)k ( 


k=l 


hr 

Ky/wTr 


= —X 


Z{x)'^af 


1=1 


iri)i + i^'^k{ipi)l 


— X 


z(i) 2aiar 


l<l<r<d 


k=2 

2 1 °o / \ / 

^(<^z)l(‘^r)l(^^) + 2 Z] ^(‘^0fe(<7>r)fc \ -^=j 


'V^r 

(3.124) 

We have obtained the expression for the asymptotic covariance (12.131) in terms of Chebyshev 
polynomials. Now we write this expression as a contour integral. Let 

Iflr 


/3:= 




make the change of coordinates x = 2yF^cos{6), y = 2y7j^cos(w), and use p2.14p to obtain that 

^ CX) 

-'^kP^{ipi)k{iPr)k 


k=l 


r'i^Fn /■2v7v / „ 

/ / </2z(a;)75r(y)rfc ( 

fc=i J-2^iJ-2^ \Wh 


7T^ 


TT^ 


dxdy 


.‘^VhJ V^V^/ \/47/ — 

/•TT ^TT 

■ / / ^/c/l^cos(/i:0)cos(A;a;)(^/(2.y/^cos0) <y9r (2Y/^cosa;) d^do;. (3.125) 

JO JO 1 


A;=l 
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Integrating by parts in 9,oj it follows that 

oo 


l^k 

—sin{k9)sin{kio) 


.k=l 


k 


^ 2 r'^ 

-'^kl5^{ipi)k{^Pr)k = ^ {2^/YlCos9) {2^cosuj) 

^ k=l ^ 

X {2y/Yisin9){2y/^sinLij)d9dijj. 

To evaluate the inhnite sum above, recall that for z G C with \z\ < 1, we have 


CO u 


\n{l-z) = - Y,— 


k=l 


(3.126) 

(3.127) 


Noting that /3 < 1, using (13.1271) . it follows that 

oo 

T 


jjk 

yy —sin{k9)sin{kijj) 


k=l 


1 ^ (5^ 
4 ^ k 

k=l 

1 - 

'4 . 

1 

"4 L 


ikO —ikO 
e — e 


ikuj —ikuj 

e — e 


- In f 1 - W In f 1 - _ in f 1 - ) + In f 1 - 


In 


(^1 -(1 -/Jeh®-'^)) -In (^1 -(l -/3eh^+^)) . (3.128) 

je 


Making the change of coordinates z = w = and recalling that (3 = this 

can be written as 

oo 


/ 3 ' 

yy — sin(/c6() sin(fc 


fc=i 


k 


'W = 


(l _/3ehe-a;)) (1 _ 

4 [ (l - /3eh®+‘-')) (l - /3eh'^+‘^)) 


--In 




1 — -^ZW ) (1 — -^ZW 
lllr J \ "/nr 


1 — '^zw ) (1 —'^zw 
lllr J V lllr 


--In 


-In 

2 


17;,. - zwp 

4 """ \.\llr - zw\^ 
Hr - ZW 


Hr - ZW 

Combining (13.1261) . (I3.129p . and noting that 

(2y^sin6() (2Y/^sinti;) d9duj = ^1 — ^1 — dzdw, 

it follows that 

^ OO 


(3.129) 


k=l 

2 

vr 




Hr - ZW 


= H \'W\'^ = Ir 

Jmz > 0 Cm re > 0 


Hr - ZW 


H 


Ir 


1 - Ch 1 - ^ ) dzdw. 
Z^ J \ yjZ 


(3.130) 
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Compare (13.1241) to (13.Sp . Using (I3.130p . (13.131) . (I3.14p and (13.91) . it follows that the covariance 
can be written as 


lim Cov{AA^^^ [(fr]} 


. Hr 




1 OO ^ X i 


2 

TT 


‘p'i{z + — )^r{w + —) i-log 


lr\ 1 


w / 2 tt 


Hr - ZW 


\z\^=H \w\^=-ir 


Hr - ZW 


1-^ 


Jm z > 0 Urn rc > 0 

7Zr(cr2 _ 2) r2^ \ipi{X) 


x (i _ 2^) dwd. + aK ~ r"' -As^Ldx -0!:Ld^. 

^ ' ^y^^Hlr J-2^ \lAji - A2 7-27> i/ A-fr “ 


3.2 The Bilinear Form 

The main goal of this section is to prove lemma [221 to which we now turn our attention. Begin 
with the following definition. 

Definition 3.5. Let M be a Wigner matrix satisfying (II. ip . and let pQ-P be the projection 
matrices defined in (12.61) and p2.10l) . For polynomial functions /, g( : M —)• M, define 


(/> 9) 


Ir.n • — 




j,k£BinBr 


= -E 
n 


Tr 


L- • 5(mW)pW| 


(3.131) 


The large n limit of (/, g)ir,n exists for polynomial functions because all moments of the matrix 
entries of M are finite. Then lim^^oo(/) 9)ir,n = {f,g)ir, where {■,-)ir is the bilinear form defined 
in definition 12.31 

We compute the bilinear form {f,g)ir of definition 12 . 31 for monomial functions f{x) = x^, g{x) = 
x^. We will also consider the random variables n~^TT{P^^^ and prove 

their convergence almost surely to the non-random limit described in lemma (231 Some results and 
techniques from free probability theory will be used. See [T] for the relevant background concerning 
noncommutative probability spaces, asymptotic freeness of Wigner matrices and for the definition 
and properties of the multilinear free cumulant functionals Hp, for p > 1. The notation is chosen 
to agree with that text. 

Regard the matrices M,P^^\P^^'^ as noncommutative random variables in the noncommutative 
probability spaces (Mat„(C), Ef^Tr]) and also (Main(C), ^Tr{-}). Since M is a Wigner random 
matrix and {P^^\ are deterministic Hermitian matrices, it follows from part (i) of theorem 

5.4.5 in [T] that M is asymptotically free from {P^^\ with respect to the functional n“^ETr(-). 
Also, it follows from part (ii) of theorem 5.4.5 in [T] that M is almost surely asymptotically free 
from {P^^\ with respect to the functional n“^Tr(-). The set of all non-crossing partitions 

over a set with p letters is denoted below by NC{p). An important consequence of the asymptotic 
freeness and almost sure asymptotic freeness of these matrices is that mixed free cumulants of M 
and {P^^\ vanish in the limit, with respect to both functionals, see theorem 5.3.15 of [T]. 
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Therefore, letting denote a product of free cumulant functionals corresponding to the block 
structure of the partition vr, it follows that 


(x^ x^) 


Ir — 


lim —E 

n—)-oo Tl 


lim -E [Tr{pWMpW • • • pWMM 
n^oo n L 

K„{p^^\ M, P^^\ ■ ■ ■ , M, p('\ p(^\ M, • • • ,M,P^^^) 

neNCi2{k+q) + l) 

Y ^c.,(M)^c^,(P«,... ,p(^’^),... ,pW), 

TTi € NC{odd), 1 X 2 G NC{even) 

TTi U 7r2 G NC{2{k + 5 ) + 1) 


(3.132) 


and also that almost surely 

Inn ^Trjp^^) (^P^MP^^^yP^^^j 

= Y At^(P('\M,pW,--- ,M,PW,PW,M,--- ,M,PW) 

7reAfC(2(fc+5)+l) 

,pM). 

TTi G NC{odd), 1 X 2 G NC{even) 

TTi U 7r2 G NC{2{k + 5 ) + 1) 

(3.133) 

Above NC{odd), for example, denotes the set of non-crossing partitions on the odd integers in 
the indicated set. Since the calculation of the joint moments in each non-commutative probability 
space (Mat„(C), n“^ETr) and (Maf„(C), n“^Tr) is identical, we make no distinction between their 
free cumulants. Lets denote by NCP{p) the set of all non-crossing partitions over p letters which 
are also pair partitions. Recall that NC{p) is a poset, the notion of partition refinement induces a 
partial order on NC{p), which will be denoted by vr < cj if, with vr, cr G NC{p), each block of ix is 
contained within a block of a. Now a notion of the complement of a partition will be developed. 

Definition 3.6. With ix G NC{pi), define the non-crossing complement ix^ G NC{p 2 ) to be the 
unique non-crossing partition on p 2 letters so that vr U vr'^ G NC{pi + P 2 ), and a < ix^ for all other 
a G NC{p 2 ) satisfying ixU a G NC{pi + P 2 )- 

Since the limiting spectral distribution of M is Wigner semicircle law with respect to the 
functional n“^ETr, and almost surely the Wigner semicircle law with respect to the functional 
n“^Tr, we have that K 2 {M) = 1 and Kp{M) = 0 for p ^ 2. It follows now that 

{x^, x'^)ir = 0, if k + q is odd, (3.134) 


and also that almost surely 

Im -Tr |(^P(^)MPW)^ (^P^MP^yj = 0, if k + q is odd. 


(3.135) 
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Supposing then that k + q is even, and continuing the calculation, 


■K 2 &NCPleven) g NC{odd) 

TTi U 7r2 G NC{2{k + 5 ) + 1) 

= E E -c„(p'",-",p‘''>) 

^ 2 &NCP{k+q) g + q+l) 

TTi < vrj 

E n 

■K 2 &NCP{k+q)i=l PdieSi 


(p(0,... 


(3.136) 


where ttJ = {S'!, • • • , 5'|7rc|} are the blocks of the non-crossing complement of a given partition. We 
have used the complement partitions to write the sum of the free cumulants over the partitions of 
the projection matrices into a product of joint moments of the projection matrices. 

Similarly, with respect to the functional n“^Tr, we have that almost surely 

= E E ,PW) 

n 2 &NCP{even) g NC{odd) 

VTi U 7r2 G NC{2{k + (/) + !) 

= E E .pM) 

n2eNCPik+q) g + g + 1) 

TTl < vrf 

E n„>i“ n p. 137 ) 

■iT 2 &NCP(k+q)i=l Pii'l&Si 


Recall that the non-crossing pair partitions are in bijection with Dyck paths, NCP{k + <?)—?• 
P(fc+g). Thus the computation for each functional reduces to counting Dyck paths. The number of 
Dyck paths (/i(0), • • • , h{k + q)) with h{k) = j is 


k 

k+j 


k 

k+j+2 


Q 

s±i 

2 


q 

i+j+‘^ 

2 


(j + 1 )^ fk + ^\f 9 + 

(fc + i)(g + i) V^y 


Note that lim^^oo ?^~^Tr for any a,b > 1. Also note that below the 

partition 7rf depends on the Dyck path d G (which corresponds to some non-crossing pair 

partition). Also note that by Ivrj^l we denote the number of blocks of Suppose for now that both 
k, q are even integers. 


The height of the path at h{k) must be even, say h{k) = 2j. Those blocks which consist only 
of the matrices P^^ will contribute a factor of to the product of joint moments. The number of 
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blocks which contain only the matrices corresponds to the number of down edges of the path 
in the first k steps. Denote by u the number of up edges and d the number of down edges of the 
path up to step k. Then u + d = k and u — d = 2j, which implies that d = k/2 — j. The number 
of blocks which contain only the matrices is equal to the number of up edges of the path in 
the final q steps. This number corresponds to the exponent on the factor 7 ,. in the product of joint 
moments. Denote now by u the number of up edges and d the number of down edges of the path 
in the final q steps. The u + d = q and d — u = 2j, which implies that u = ( 7/2 — j. The remaining 
blocks of the partition contain projection matrices of mixed type and will contribute a factor jir 
to the product of joint moments. Since the total number of blocks in the partition is + 1, the 
number of factors of 'jir in the product of joint moments is 2j + 1. Partitioning the Dyck paths 
into equivalence classes based on the height h{k), we get that 






Ir — 


E n„i- n 

rfeD(fc+,)i=i p0)e5i 


E E 

^=0 dGD(fc+,) 


2-i 2 -^ 2i+i 
7 r hr 


h{k) = 2j 

( 2 j + 1 )^ f k + l\ f g + l\ 2 j+i 


E (fc -g l)(g + 1) j (^ g+2j+2 JT; 7r hr > 


j=0 


and also, almost surely, 


lim — Tr 

n—)-oo Tl 




= E np“ n 

rfeP(fc+,)j=i P7)e5i 

2 

= E 

2~J 2j+i 

j=0 



II 

to 

k 

2 

( 2 j + 1)2 f k + l\f q + l 

- E 

j=0 

{k + i)(g + 1 ) 


2~^ 2 J 2 j + l 

h 7 r hP ■ 


Now suppose that both k, q are odd. The height of the path at h{k) must be odd, say h{k) = 
2j + 1. Similar to the even case, the number of blocks which consist only of the matrices P^^'^ 
equals the exponent of h in the product of joint moments. The number of blocks which contain 
only the matrices corresponds to the number of down edges of the path in the first k steps. 
Denote by u the number of up edges and d the number of down edges of the path up to step k. 
Then u + d = k and u — d = 2j + 1, which implies that d = {k — l)/2 — j. The number of blocks 
which contain only the matrices P^^^ is equal to the number of up edges of the path in the final q 
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steps. This number corresponds to the exponent on the factor 7 ^ in the product of joint moments. 
Denote now by u the number of up edges and d the number of down edges of the path in the final 
q steps. The u + d = q and d — u = 2j + 1, which implies that u = {q — l)/2 — j. The remaining 
blocks of the partition contain projection matrices of mixed type and will contribute a factor of jir 
to the product of joint moments. Since the total number of blocks in the partition is + 1, the 
number of factors of 7 ;^ in the product of joint moments is 2j + 2. Partitioning the Dyck paths 
into equivalence classes based on the height h{k), we get that 




1 , 


dG-Df^k+g) ^=1 P^^'^GSi 


k-l 

2 


E E 

deD 


2j+2 

h 7r " h; 


(fe+g) 

h(A:) = 2j + 1 


h-l 

2 


(2j + 2) ( g + 1 ^ 2j+2 


i=o 


and also, almost surely, 

lim — Tr 

n—)-oo Tl 




kjl 


E ^Tr{ n 

deD(fc+g)i=l pO')eSi 


fc-1 

2 


E E 

^=0 d e 

h{k) = 2j + 1 


2j+2 

h 7;/ 


k-l 

2 


(2j + 2)^ ( k + 1^ ( q + 1 ^ 2j+2 

= E + + 7 . T. . 

Now for polynomials f{x) = 5'(^) ~ YlT=o^j^^^ have by linearity that 


p m 


{f,g)lr = '^'^aibj{x\x^)lr. 
i=0 j=0 


(3.138) 


The intersection of countably many events, each with probability 1, occurs with probability 1. 
There are only countably many polynomials with rational coefficients, so we have proved that the 
random variables 

-Tr{p(') 

n 

converge almost surely to the same, non-random limit given by the right hand side of (13.1381) . 
whenever /, g are polynomials with rational coefficients. 
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In the next proposition the bilinear form {f,g)ir is diagonalized. 

Proposition 3.7. The two families and {Uq"'}’^Q of rescaled Chebyshev polynomials of 

the second kind are biorthogonal with respect to the bilinear form 113.131]) . More precisely, 


= 4 , 


Proof. Since {x^, x'^)ir = 0 if A: + g is odd, it follows by linearity that 

{U^\U^^)ir = 0, if A: + O' is odd. 

We begin by computing U]g)ir and U]q^^)ir. We obtain 






y/Tl 

7,=n \ P J 


2(?-2p+l [k + j + l)[ 


P 


(3.139) 


(3.140) 


p=0 

7; ^7r- ^ (~l)^7f(2j + 1)^ I' 2k + 1 ^ ^2(7 — ^ 20^ — 2p + 1 ^ ^k-j^^q-p-j^2j+i 


P + J + lP 


7h 


1 


2A: + 1 


D 2 J' + 1)^(7 


7 2A: + 1 \ 

^ (_l)P(2g-p)! 

\k + j + 17 

^ p!(g - p + J + l)!(g - p - 7)! 


—J — i ‘^3 1 
7/ ' It ^ hr 


(3.141) 


and 


(( 


X 


\ 2 /cH-l TTlr 


2y/7l' 


Ufg+l)lr 




y p'? P + A ,2k+l^^2q-2p+U 

V3i tr'n \ P J 


p=0 

"o 2 ^ ^ 3 

h 7r 


(—1)^7^ (2j + 2)^ / 2A; + 2 \ /2(7 —p+l\/ ‘2g — 2p + 2 
2^ + 2 ^ ^ 2q-2p + 2 VA: + j + 2^ V P 7 U “ P + i + ' 


\ —7 —Tl— 


7 V? - p + 7 + 2/ ' 


P,-p-j^3j+‘^ 


hi 


-1 -1 p 

22 


7; ^7r 
2A; + 2 


■E(27 + 2) 


i=o 


2 / 2A: + 2 \ 

U + j + 27 


(l)P(2g p+l)! 


V- 

^ p!(g - p + 7 + 2)!(g - P - j)! 


-1 -i 2j+2 
7i Ar ^hr ■ 


(3.142) 
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Denote by 


Q-J 


{-l)P{2q-p)\ 


^ p'-{q -p + j + i)!(9 -p- j)'-'' 
g-j 


H2{qJ) = 


p=0 


{-l)P{2q-p+iy. 

pKq - p + j + “^V-iq - p - j)' 


Then 




-1 _i p 

9 oft/ 


i=o 


' 2 Vt^ 

It follows from (13.143113.144jl that 

{2q)\ 


H2{q,jhi 


Hi{q,j) = 


{q - jy.{q + j + 1)'. 


= 2 F 1 


-{q- j),-iq + j + ^) 

-2q 


;ih 


(3.143) 

(3.144) 

(3.145) 

(3.146) 

(3.147) 


H2{q,j) 


(2g + l)! 

iq-jy{q+j + iy 


-(9-j)>-(9 + J + 2) 

-i2q+l) ’V’ 


(3.148) 


where 2-^1 is a hypergeometric function. See [3] for the definition of hypergeometric functions. 
Below let {x)n = x{x + 1) • • • (x + n — 1) denote the rising factorial. By the Chu-Vandermonde 
identity (see e.g. [3]), it follows that 


Hi{q,j) = 
H2{q,j) = 

Therefore, for A; = 0,1, • • • 
0. With k = q we obtain 


{2qy 

{3-q + l)g-i _ J 

f 0 

0<j<q 

{q-jy{q + j + 1 )! 

{~‘^q)q-j 1 

1 2q+l 

J = q 

(2<? + l)! 

{j-q + i)g-i J 

f 0 

0<j<q 

{q - jy{q + j + F)'. 

{-2q-l),_, ] 

1 2q+2 

j = q 

,q 1 , we get that ,U(()ir - 

0 and also (( 2 ^) 


(3.149) 

(3.150) 


'2q+l/ 


1 


2k+ 1 


k 

^(2j + lf( ^ \n. (h 

7 —n 


3=0 

{2k + lf (2k + I 
2k + I 

ihr 




(3.151) 
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and 


// ^ \2k+l rpr \ 


-1 _i i. 

2 2 ^ 


E( 2 i + 2 )' 


7; 7r 

2A: + 2 ^ 
i=o 

_i -U2k + 1)^ f2k + l 


2k+ 2 
k + j+ 2 


H2{k,jhi ^7r 


-io,-io, 2 i +2 


= 7; ^7r 

2fc+2 
Hr 


2k+ 1 \2k + 




-2A;+1 • 


Thus, for k < q, 
and for k = q 


yfWr 

{Ul, U^^)ir = 0, (C/^V, C/2";+i)z. = 0, 


2 A :+1 


(f^ 2 fc> ^Tk)lr - ’ 


(^2fc+l’^2fc+l)*^ ~ ((o ^2fc+l)^»’ ~ 


7z. 


2fc+2 


2^i' ’- 2 fc+i/- ^ 2 fc+r 


(3.152) 

(3.153) 

(3.154) 

(3.155) 


This completes the proof of proposition 13.71 which is the diagonalization part of lemma [2l5l 

Remark 3.8. Previously we have shown that whenever /, g are polynomials with rational coefficients, 
almost surely (a.s.) 

hm iTr |P«/(M«) • 5 (mM)pM| = {f,g)ir. 

n—>oo n y. > 

The Chebyshev polynomials have rational coefficients, so it follows from the above argument 
that a.s. 


1 


v'T7r n 




Hr Nfc +1 


^fiar 


(3.156) 

□ 


Now the bilinear form {■,-)ir will be extended to functions other than polynomials. For this 
part of the argument, the bound on the variance of linear eigenvalue statistics in 13.31 is essential. 

Proposition 3.9. Let f,g & Tig for some s > |, i-e. for some e > 0, 

/ oo poo 

\f{t)\Hl + \t\f^"dt<oo, / |?(t)|2(l + |t|)3+^dt <oo. (3.157) 

-OO j —OO 

Then the limit of {f,g)ir,n (see definition ] 3. 5\} as n ^ oo exists and 

r2v77 /■27> 


{f,9)lr = 




f{x)g{y)Fir {x,y)y/A'ji -x^\/ - y^dydx, (3.158) 
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and also, almost surely, 




1 r‘iVTi - - 

;- / / f{x)giy)Fir{x,y)\/A'yi-x‘^y'A'yr-y‘^dydx, 

"lllr J—2^PTi J —2^Py^ 


Fn’^iar J-2^J-2^ 


(3.159) 


where the kernel Fir{x,y) is given by id.8ti\) . 

Proof. First it will be argued by approximation that (•, ■)ir can be extended to the class of functions 
^3 ,,, and then the bilinear form will be explicitly computed. It will be sufficient to approximate 

rrr-rn 

/, g below by truncated polynomials with rational coefficients in "Hs , because of the estimate 13.31 . 
Recall that functions of the Schwartz class are dense in T-Lg, so after a triangle inequality argument 
it is in fact sufficient to suppose that f,g ^ 5(]R). Let h S be a function so that h{x) = 1 
for X € [—3,3], h{x) = 0 for x ^ [—4,4] and is smoothly interpolated in between. Note that with 
overwhelming probability, the eigenvalues of the submatrices concentrate in the support of g,sc- As 
a consequence we may suppose that f,g are supported in [—3, 3]. We give a density argument. 
It is sufficient to argue that \\hf — hpj\\3_^^ and \\hg — hqj\\3_^_^ converge to 0 as j —)• oo, where 
{pj},{qj} are appropriately chosen sequences of polynomials with rational coefficients. Note that 
hf = f and hg = g. We now focus on estimating ]]/ — hpj\\3_^_^. Since / is a Schwartz function, we 
have that / £ 7i2. We note that 

/ CXD POO 

\m\"{l + \t\)^+^dt< / \f{t)f{l + \t\)Ut, (3.160) 

-oo J —OO 


so it will be sufficient to approximate / in the larger jj • jj 2 norm. Also, since 


/ OO ^ 

l/(t)l^ (I + [tD^dt < Const 

-OO 


\f{t)\‘^dt + 


/ OO 

t^\M 

-oo 


dt 


(3.I6I) 


we only need to approximate the two terms on the right hand side. Consider polynomials {pj} 
with rational coefficients so that sup_ 4 < 3,<4 j/"(x) — pj{x)\ —^ 0 as j —>• oo. Then denote by 
Pj{x) = j^^Pj{t)dt, and Pj{x) = j^^Pj{t)dt. As a consequence of Parseval’s theorem, it will be 
sufficient to show that 


11/- %IIl 2 ([_ 4 , 4 ])-S'0 and ]]/" - (/ipj)"l]i 2 ([_ 4 ^ 4 ])0, as joo. (3.162) 

But observe that 

II/" - (%)"IIl 2 ([- 4 , 4 ]) < II/" - ^PjIIl 2([-4,4]) + + 2/ip/|i^2([_4^4]). (3.163) 

The first term on the right hand side converges to 0 because of the uniform approximation. 
Noting that h'{x) = 0 and h"{x) = 0 on (—3,3), and also that pj and pj converge to 0 uniformly 
on [—4, —3) U (3,4], it follows that the second term on the right hand side converges to 0 as well. 
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Finally we observe that 


II/- 


([-4,4]) 


< J h^{x) 


j Jfi^) - h{x)pj{x)fdx 

J J [f"{u) - pj{u)]dudt 

< Const- ( sup \f''{u)—pj{u)^ 
V-4<u<4 V 


dx 


(3.164) 


It follows that 11/ — h/jl|| 2([_4 4 ]) —> 0 because of the uniform approximation. This completes the 
approximation argument, so we now turn toward computing the bilinear form. 

Setting 

/fc = X- / f{x)U^‘{x)^/4'yl - x^dx, gu = - - / g{y)Ul'' {y)^/- v'^dy, (3.165) 

2vr7/ J-2^i J-2^ 

it follows that 
(/, g)ir 

OO OO 

= ^gpU:^^{x))ir 

k=0 V=0 

OO OO 

k=0 p=0 

OO ^k+1 

^ ^ -^kffk ^/2 f^/2 
k=0 h 7r 


1 


C\Cr 


f{x)g{y) 

-2./^ 


It also follows, using p.l56p . that a.s. 




7« 


fe+i 


Lfc=0 


k/2 k/2 

h Ir J 


V47; -X^s/47 r - y'^dydx. 


CVTl /■2v7v ryf 

4 2 / / fix)g{y) 

Air^-far J-2^i J-2^ zcL -iV 


k+l 


.k=0 


k/2 k/2 
h Ir 


\/ 47; — x^\/ 47r — y^dydx. 

(3.166) 


Proposition 13.91 follows. □ 

Lemma [2.51 follows now from propositions 13.71 and 13.91 This also completes the proof of theorem 
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4 Proof of Theorem 12.2 


It is enough to show the case of d = 2, i.e. the limiting covariance of and Then 

by symmetry property, one will get the limiting covariance of and N'n'^° [fp], 1 < I < p < d. 

Let U{t) ,U{t),Un{t),Un{t) heU^^\t), 17^^^ (t ), Un^ (t), Un^ (t) defined in (l,'I.16ll,'I.17l) respectively. U (t) 
and U{t) are unitrary matrices and 


U{t)U*it) = U{t)U*{t)=I, \Ujk\<l, 
By Remark 3.3 in m, we have the following bounds 




k=l 


Var{Un(t)} < C'(cr6)(l + |t|^)^, 

(4.1) 

Var{u„(t)} < C'(cr6)(l + |t|^)^, 

(4.2) 

CXD \ 2 

Var{Xi)°(t)} < Ciae) (1 + |t|^)l^i(t)| dtj , 

(4.3) 

Var{X2)o^^)| < ^ _ 

(4.4) 


Let rc be a linear combination of random variables and N'n^°[<P 2 ], and Zn{x) be the 

characteristic function of w, i.e. 


r(2)or 


So 


w = Zn{x) = E{e*"“'}. 

Zn{x) = 1 + r ZUt)dt; Z^x) = iE{we^^^}, 

Jo 

By Cauchy-Schwarz inequality and (|4.3ll4.4p 

/ CX> 

(1 + |t|^)|^i(t)| dt, 

-OO 

Fourier inversion formula /(A) = f dt givies us that 

/ OO poo 

^(tX(t)dt, A/'XX] = / 

-OO j —OO 


(4.5) 

(4.6) 

(4.7) 


So 


where 


w= a^{t)Un{t) + I3(p2{t)u°n{t)dt, 

J —OO 

/ OO poo ^ 

^{t)Yn{x,t)dt + i/3 i^{t)Yn{x,t)dt, 

■OO J —OO 

Yn(x, t) = E{un{t)en{x)}, Yn{x, t) = E{u° (t)en(a;)}, en{x) = e"' 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


36 









By Cauchy-Schwarz inequality, 


F„(x.t)l < E{K(*)I) < + |i|“), 

|F„(i,t)| < E{|<(t)|} < C'/VeXl + 1*1’). 


and 


I — 

dx 


Yn{x,t)\ 


= |E{a<7VW°[(^i]e„(x) + /3<A/'P°[(/?2]en(a;)}| 

/ OO 

(1 + |t|^)(|a^i(t)| + \l3(p2{t)\)dt. 

-CO 


And 


O 

—Ynix,t) =E{u’^{t)e°^{x)} = ^ 


E{Wjk<^n] 


where 

= Ujk{t)e°^{x). 

Recall Djk = djdMjk, /3jk = (1 + djk)~^, 


DjkUab{t') — ^j,k£Bi'i‘(djk{Uaj * “ 1 “ ^bj * Uak(t')]^ 

d^jk^abid) ^j,k&B2^^jk{Uaj * U^kit) Y Ujjj * C^afe(^)]' 


and 


(4.12) 

(4.13) 


(4.14) 


(4.15) 


(4.16) 

(4.17) 


DjkSn (x) 


2i/3jkxen(x)(lj^keBia(Yi)jk(^i) + li,fceiJ2/3(7^2)jfc(^2)) 

( ^oo roo ^ 

lj,keBi tUjkit)a(^{t)dt + lj^keB2 tUjkld^{t)dt 

J— OO j—OO 


(4.18) 

(4.19) 


Lemma 4.1. Assume <^ 1 , 7^2 have fourth bounded derivatives. Then 

\Djk{Ujk{t)e'f{x))\ < Ci{x,t), 0<l <5 (4.20) 

where Ci{x,t) is a degree I polynomial of \x\, |t| with positive coefficients. 


Proof. Prom (|4.16p and (I4.17h . we have 



\D]kUab{t)\,\D]kUab{t)\ < Consti\t\\ 0 < Z < 5. 

(4.21) 

(14.19D imolies 

\D^jj^en{x)\ < Const'i{l + x|^) 0 <l < 5. 

(4.22) 

These two complete the proof of Lemma 14.11 

□ 

AddIv (jA.iP with n - 

= 2 to obtain 


|F„(x,t) = 

^ + l)6,kMD,k^n} + 0(1) 

j,keBi 


= 

- (1 + djk)E{Djk<^n} + ^ E{Dj,<^n} + 0(1). 

” jAsBi ^ iesi 

(4.23) 
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where the error term is bounded by C'i{x,t) as n —>■ oo. The first term in (|4.23l) is 


t 


1 


rt 


1 , 


rt 


- Yn{t,x) -/ E{un{t-ti)}Ynix,ti)dti -E{ / Un{ti)u°^{t - ti)dtie°^{x)} 

n n Jo n Jo 

2i 

-E{xe„(x) 

n 


hunit + ti)a^{ti)dti + / tiTrP^^’‘^^Un{t)P^^’‘^^Un{ti)l5(p2{ti)dti 


The first term and the second term are bounded because of (j4.12p . The last term is bounded by 
XI(|q^II^i(^i)I + |/5||^2(ti)|)hti- And the third term is bounded by 2\t\C^/‘^{a%){l + 

The second term in (|4.23l) is 


2-a^ 


n 


+ 


ix{2 — cr^) 


TPt. Uo 


' - h)dheUx)] + 


n 


■ E|en(x) f 

jeBi '' 


tiUjj{t)Ujj{ti)a^{ti)dti 


n 


Y E|e„(x) f tiUjj{t)Ujj{ti)l3i^{ti)dp 


jeBinB2 


The first term is bounded by 2|2—cr^ | 
|/3||^2(fi)|)dfi- So 


and the second term is bounded by 2|x||2—|t| (|a||^i(ti)| + 




< C^ixp). 


Symmetricly, Yn{x,t) has similar bounds. So we conclude that the sequences {Yn},{Yn} are 
bounded and equicontinuous on any finite subset of We will prove now that any uniformly 
converging subsequence of {Yn}{{Yn}) has same limit Y{Y). 

We deal with Yn first, and by the symmetric property, we can hnd Yn- By identity, 


Unit)=ni + i Y MjkUjk{ti)dti, 

■^0 7,fees. 

(4.24) 

so 

/.V- 


Yn{x,t) = ^ / V E{WjkUjk{ti)e°^{x)}dh. 

(4.25) 

Bv aonlvine decouoline formula (lA.lD with » = 3 to (14.251). we have 


j,keBi 

+£3,,* 

J=0 ^ 

dti, (4.26) 

where 



^IJk ~ 

= 0,K2,jk = 1 +Sjk{cr'^ - 1), 

(4.27) 

^3,jk — /^3? kvA,jk — ^4? J 7^ 

(4.28) 

and kisjj, kiAjj are uniformly bounded. 

i.e. there exist constants (T 3 ,(T 4 such that 

|^3,Xxl — ^3? 1^4,jxI — ^4? 

(4.29) 

and 



ksjfcl < n ^C'3E{|fTjfcX}sup|T)|fc4>n(x)| < n "^C^ixp). 

tew. 

(4.30) 
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Let 


„(Z+l)/2 


E 


'0 — I'- 


^^+^dk^{Di^^Ujk{ti)e°nimdti,l = 1,2,3, 


f — 

<-'n. — 




(4.31) 

(4.32) 


j,keBi 


Then 

and 


Yn{x, t) — Ti + T2 + T3 + £n, 


\£n\ < ) C^{x, t) ^0, as n ^ 00 . 


Note: if Wj^s are Gaussian, then Yn{x,t) =Ti. So Ti coincide with the in Theorem 12.11 
Let 

Vn{t) = n~^E{un{t)}, Vn(t) = n“^E{u„(t)}. 

Then 


Yn 

where 


{x,t) + 2f dti f Vn{tl-t 2 )Yn{x,t 2 )dt 2 =xZn{x)An{t)+rn{x,t)+T 2 + Ts+£n, (4.33) 
Jo Jo 


A„,{t) = -— [ E{TrU{ti)Piip[{M,)Pi}dti [ E{TrU{ti)P 2 ^' 2 iM 2 )P 2 }dti, (4.34) 
n Jq b Jq 

and rn{x,t) —)• 0 on any bounded subset of {(3:,t) : x G M, t > 0}. 

Let A{t) = lim^^oo ^n(i)- From the proof of Theorem 12.11 A{t) coincides with the one in 
Gaussian case. 

Proposition 4.2. r 2 —>• 0 on any bounded subset of {{x,t) : x G M, t > 0}. 

Proof. The second derivative (1=2) is 
D%{Ujk{ti)el{x)) =/3|fcX 

{-{Wjj * Ujk * Ukk + “JUjk * Ujk * [/jfc)(ti)e°(x) 


-^i{Ujj * Ukk + Ujk * Ujk){ti)xen{x) 


tUjk{t)w^i{t)dt+ lj^k&B2 / tUjkl3^{t)dt 


+4[/jfc(ti)x^e„(x) 

-2iUjk{ti)xen{x) 


/ OO 1*00 

tUjk{t)ai^{t)dt + lj^keB2 / tUjk(J^{t)dt 

■OO J—OO 


P^j,keB2 


/ CO ^ 

m 

-OO 


' —OO 

roo 


/ OO 

^(Ujj * Ukk T Ujk * Ujk^{tJ)(y.ip\{t)dt 

-CXD 


ii * Ukk + Ujk * Ujk)it)/3^{t)dt 


}• 
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Let 


^21 = * Ujk * Ukk + 2U,k * U,k * Ujk){ti)e°^ix) 

* Ukk + Ujk * Ujk){ti)xen{x) J t2Ujk{t2)a^{t2)dt2 
+4/3jkUjk{ti)x‘^en{x){ J t2Ujkit2)a^{t2)dt2f 

-2i/3jkUjk{ti)xen{x) / t2{Ujj * Ukk + Ujk * Ujk){t2)f3(^{t2)dt2 > dti, 


T 22 = / E < ^ 4l3‘jkUjk{tl)x‘^en{x){ / t 2 Ujk{t 2 )P^{t 2 )dt 2 f 

^ {j,keBinB2 

+SI3^jkUjk{ti)x^en{x) j t2Ujk{t2)a^{t2)dt2 J t2,Ujk{t2,)l3^{t‘i)dt^ 

-2if3jkUjk{ti)xenix) / t2[Ujj * Ukk + Ujk * Ujk]{t2)/3i^(t2)dt2 > dti, 


i 

^23 = / Y1 ''^3,jj^{D‘jj{Ujj{ti)e°^{x))}dti. 

-^0 jesi 

Then T 2 = T 21 + T 22 + T 23 . In [TO], it has already shown that |T 2 i| < |t|C' 2 (x,t)ni/n ^/2 on 
any bounded subset of {{x,t) : x G M, t > 0}. Also, by Prooersition 14.11 and (14.291) . one has 

17231 < |t|C'2(x,t)ni/n3/2_ 

In T 22 , there are three types of sum, 

51 = n-3/2 ^ 

j,keBinB2 

5 2 = n-3/2 ^ Ujk{ti)Ujk{t2)Ujk{t3), 

j,k&BinB2 

53 = n-^/^ E Ujk{h)Ujj{t2)Ukk{t3). 

j,keBinB2 


By the Schwarz inequality that 

|5i|<n-3/2 ^ \u.,^t2)Ujk{t3)\<^^ 
j,keB2 

|52| < n-3/2 \U,k{ti)U,k{t2)\ < ^ 


Write 

*^3 = ^{Pl2U{h)P^2V{t2),V{t3)), 
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where 


V{t) — ni 2 ^ (^jji^))j&BinB 2 - 

||^(i)|| < 1) ||-Pi 2 ^^(i)-Pi 2 || < 1) so we conclude that S 3 < hence T 22 < |t|/n^/^. This completes 
the proof of preposition 14.21 □ 

Proposition 4.3. 

T3 = T31 + T32 + Rsix, t), 

where 


T 31 = 


IK 4 

rR 


IL E |c/jj * Ukk{ti)xen{x) J t2Ujj * Ukk{t2)a(^{t2)dt2'^ dti, 


j,keBi 


T 32 = 


ZK4 


/' E 

*^0 A U.<~\ 


¥.\Ujj *Ukk{ti)xen{x) / t 2 Ujj *Ukkit 2 )l 3 ^it 2 )dt 2 > dti. 


j,kGBir]B2 

and R 3 {x,t) 0 on any bounded subset of {{x,t) : x S M, t > 0}. 

Proof. 

iK4 


= E nD%{U,k{ti)el{x))]dt3+T3, 


j,k£Bi 


where 


^3 = ^ / E - Ki)^{D^j{Ujj{ti)e°^{x))}dti. 


j&Bi 


By Propersition 14.11 and (j4.29l) . we have IT 3 I < |t|C' 3 (x, t)ni/n^ 
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The third derivative ( 1 = 3 ) 

D%{U,k{h)el{x)) = 

* Ujk * Ujk * Ukk + 617jj * Ujj * Ukk * Ukk + ^Ujk * Ujk * Ujk * Ujk){ti)e^{x) 


+6{6Ujj * Ujk * Ukk + ‘^UjkUjk * Ujk){ti)xen{x) ( / tUjk{t)QA^{t)dt 
+lj,fceS2 / tUjkPif2it)dt 

+12i{Ujj * Ukk + Ujk * Ujk){ti)x^en{x) tUjk{t)a^{t)dt + Ij^k&B^ J tUjkl3^{t)dt 

+6{Ujj * Ukk + Ujk * Ujk)iti)xen{x) (^j t{Ujj * Ukk + Ujk * Ujk){t)a^{t)dt 

+lj,fce-B2 J t{Ujj * Ukk + Ujk * Ujk)(d^{t)dt^ 

-8Ujk{ti)x^en{x) tUjk{t)a^{t)dt + lj^keB2 j tUjk/d^{t)dt^ 

+12iUjk{ti)x'^en{x) tUjk{t)a^{t)dt + lj^k&B2 J tUjkP^it)dt 


Jm 


33 * ^kk “ 1 “ Ujk * Ujk)(i)(y-^l{i')dt lj_^kGB2 I ^{Ujj * Ukk “ 1 “ Ujk * Ujk)P^2(t)dt 


jm. 


+ 2 Ujk{ti)xen{x) j t{ 6 Ujj * Ujk * Ukk + ‘ 2 Ujk * Ujk * Ujk){t)a^{t)dt 

+^j,keB2 J K^Ujj * Ujk * Ukk + ‘^Ujk * Ujk * Ujk){t)fd^{t)dt 

So any term of 


}■ 


( 4 . 35 ) 


^/^4 

Gn^ 


[ Y 1 nD%{Ujk{h)eUx))}dh 
Jo 


j,kGBi 


containing at least one off-diagonal entry Ujk or Ujk is bounded by ( 73 ( 3 :, t)ni/n^. Let Rz{x,t) be 
the sum of T 3 and these terms. Then |ii 3 (a;, t)| < C' 3 (x, t)ni/n^-|-|t|C' 3 (x, t)ni/n^. So two terms in 
(I 4 . 35 P containg diagonal entries of U and U only left contribute to T 3 . They are T 31 and T 32 . □ 


Let 

r-2v^ 
'-2^71 

By Wigner semicircle law, one has 


I /■2v^ . ^ ,_ 1 ^_ 

v{t) = - - / v{t) = - - / 1/472 - X^dX. 

27r7i 7-2,77 2^72 J-2^ 


lim Vn{t) = livit), lim u„(t) = 72 u(t). 
n—)-oo n—)-oo 


Then 


{v * v){t) = --—2 / e*‘^/rV47i - 

2vr7i 7-2771 




TTtjf 7-2771 1 / 471 - 72 


d^i, 


( 4 . 36 ) 
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(4.37) 


j _ 1 r 

(v * v)(t) = --—2 / e*‘^/iA/472 - / 

2vr7| J-2^ vrt7^ J_ 


1 Jt,j72- 




Let 


Denote 


Proposition 4.4. 


7rt7| J-2^ ^472 - 

l(t)= [ {v *v){ti)dti, I{t) = [ {v*^{ti)dti. 

Jo Jo 

1 /■ 2 \/^ 27 / - 

B^i = —2 / = 1,2. 

71 - 7 / 4-27^ V^7; - 

T 31 iK 4 xZ{x)I{t)a'ylB^j^, 


dfi. 


T 32 iKAxZ{x)I{t)l3'yl2B^^. 

uniformly on any bounded subset of {{x,t) : x € M, t > 0}. 
Proof. The proof of (I4.40p can be found in [T^. And 


(4.38) 

(4.39) 


(4.40) 

(4.41) 


T 32 = 


IXK4 


t 2 ^{Ujj{ts)Ukk{tl — tz)Ujj{t/fjUkk{t 2 — ti)en{x)} 




j,keBinB2 

X I3^{t2)dtidt2dtzdti 
rt rti r ft2 


n ci r rt 2 

/ / t2E{Vn{t3,t4)Vn{tl - ta, ^2 “ t4)e'^{x)} /3l^{t2)dt4dt2dt3dti 

J Jo 

rt rti r |•t 2 

+iXKiZn{x) / / / t2^{Vn{t3,t^)Vn{tl-t2„t2-ti)}l5L^{t2)dtidt2dt3dti 

Jo Jo J Jo 


where 

Then 
and 
where 

Proposition 4.5. 


Vn{tlM)=n^ ^ Ujj{tl)Ujj(t 2 ). 
ieBinB2 


(4.42) 

(4.43) 


\&{Vniti,t2)Vn{H,t4)el{x)]\ < 4E{|x° (ti, t2)|} + 4E{|(ts, t4)|}, (4.44) 

E{'^;n(^l,^2)^’n(^3,^4)} = hn(t 1, t2)hn(ts,^4) + E{Vnitl,t2)v°{t3,U)}, (4.45) 

Vn{tl,t2) ='\^{Vn{tl,t2)}. (4.46) 


Vn{tl,t 2 ) = 'yi 2 v{tl)v{t 2 ) + o(l), 
uniformly on any compact set o/M^. 

Proof. Indeed, K{Ujj{ti)Ujj{t 2 )} = v{ti)v(t 2 ) + o(l) uniformly in 1 < j < n and ti,t 2 from a 
compact set of , which follows from 


EUjj{t) = v{t) + o(l), Var{Ujj{t)} = o(l), 
IKUjjit) = v{t) + o(l), Var{Ujj{t)} = o(l) 


(see e.g. [H]). 


□ 
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So the limit of T 32 is 


pt poo 

ixK^Z{x )^\2 / v*v{ti)dti / t 2 ( 3 i^{t 2 )v *v{t 2 )dt 2 = ixKiZ{x)'yl 2 ^i'f^)PB^^. 

Jo j — CO 


□ 


So if Y{x,t) = lim^^oo i)j then Y{x,t) satisfies 
rti 

1—00 Jo 

Therefore, if let Y*(x,t) be the solution of 

rti 


/ oo pt\ 

dti / v{ti-t 2 )Y{x,t 2 )dt 2 = xZ{x)[A{t)+iKAl{t){a'ylB^^+j3^‘l2B<f2)\ 

-OO Jo 


then 


/ OO pt\ 

dti / v{ti - t 2 )Y{x,t 2 )dt 2 = xZ{x)A{t), 

-OO J 0 


iK4xZ{x) r 2 d fl 2 d 1 e®^'^( 27 i — A^) 


Y{x, t) = Y* {x, t) + ""'"^ 2 '"' [«7i+ ^7 i2^¥’2] [ /--^ 

27r7i V47 i - A2 

Symmetrically, 


dX. 


Y{x, t) = Y*{x, t) + 

27r72 


e*‘^(272 - A2) 

/-2V7J 472 - A 2 


hA. 


(4.47) 


(4.48) 


Therefore 
Z 


/ OO ^00 ^ 

^(t)y(x, t)dt + i ;8 / ^{t)Y{x,t)dt 

-OO J —00 


= —a;l/Z(x) — a 


K4xZ{x) 

2-K-il 


^(0 +/ 3712 -BVP 2 ] [ 

I J — 


;?S(t) [«7?2 ^v^i +/37i^v^2] 

^'^'12 J— 00 


2YY e*‘^( 27 i - A2) 
-2^71 V 471 - A2 
/■2V77 e*^^(272 -A2) 

I - 

'-2^7^ V 472 - A2 


dXdt 


-dXdt 


= -xVZ{x) - Q, 2 ^yb i 72 ^$i - af 3 xZ{x)'yl 2 Bip-,B^^ - ^^^ 7172 ^ 2 ^ 


= —xVZ{x) — XK 4 Z{x) 


VI 

«^7?^ + a/djf^B^iB^^ + /3^7|^ 


(4.49) 


where 


y = a2Var(Gi) + 2a/3Cov(Gi, G 2 ) + p^Var{G 2 ), 

and Gi, G 2 are the random variables in Theorem 12.II with d = 2. 
Therefore, 

lip Cov(AAW°[v9i],Arp)°[(^2]) = 


Cov(Gi,G2) + „Z9^?2 f V’iilJ)^^==dn 




2vr27i72 J-2^ a/471 - /r2 7-2^77 vT2^T^ 


zd/i. (4.50) 
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By symmetry, for any 1 < / < p < n, 


Cov{Gi,Gp) 


Cov{Gi,Gp) + 







27/ - 


dA 



¥^p(/^) 


V^47p -/r2 


Appendices 

A Decoupling Formula 

Theorem A.l (Decoupling Formula). /? J llO^ Let ^ be a random variable such that ^ qq 

for a certain nonnegative integer p. Then for any function / : M —)• C o/ the elass (7^+^ with 
bounded derivatives f^^\l = l,...,p + 1, we have 

E{f/(f)} = E^E{/''>(0}+£p- (A.1) 

1=0 

where ki denotes the Ith cumulant of f and the remainder term Sp admits the bound 

kp| < C,E{|Cr+nsup/(^’+i)(t), C, < l + (3 + y+^ _ (A 2) 

teK [P+^)- 

If f, is a Gaussian random variable with zero mean, 

EU/(0}=EU'W/'(0}- (A.3) 
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